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Abstract 

The newly developed iterative method based on Green function defined by 
quadratures along a single trajectory is combined with the variational method to 
solve the ground state quantum wave function for central potentials. As an exam- 
ple, the method is applied to discuss the ground state solution of Yukawa potential, 
using Hulthen solution as the trial function. 



PACS: ll.10.Ef, 03.65. Ge 



1. Introduction 



Recently an iterative solution of the ground state for the one dimensional double-well 
potential is obtained[l] based on the Green function method developed in ref.[2]. This 
Green function is defined along a single trajectory, from which the ground state wave 
function in N-dimension can be expressed by quadratures along the single trajectory. 
This makes it possible to develop an iterative method for obtaining the ground state wave 
function, starting from a properly chosen trial function. The convergence of the iterative 
solution very much depends on the choice of the trial functionfl]. Therefore the key of 
this new method is to find the best trial function as the starting point of the iteration 
process. 

On the other hand the variational method is a well-known approximate approach 
to solve Schroedinger equation starting from a trial function with a set of variational 
parameters. By minimizing the expectation value of the Hamiltonian respect to the set of 
parameters of the trial function, the variational solution could be obtained. If the staring 
trial function is properly chosen the solution could be a very good approximation for the 
equation. However, it is difficult to further improve the accuracy of this approximate 
solution. 

In this paper the two methods are combined to get a veriational-iterative solution. The 
result is obviously more accurate than the variational one. It has been shown in ref.[l] 
that the final iterative solution is convergent if the trial function is properly chosen. As 
an example the variational-iterative method is applied to solve Yukawa potential, using 
the solution of a Hulthen potential as the trial function. The obtained result up to the 
second order is compared to those based on the pure variational procedure, starting from 
the same trial function. The solution of the combined method is improved significantly. 

In Section 2, a brief introduction is given about the Green function method based 
on the single trajectory quadrature and the corresponding iterative formula for central 
potentials. The variational-iterative solution for Yukawa potential is given in Section 3. 
The numerical results and discussions are given in Section 4. 

2. Green Function and Its Iterative Solution 
for a Central Potential 

1. Green function 

We discuss a particle with unit mass, moving in a central potential V(r) 
state wave function \I/(r) satisfies 

[--V 2 + V{r)]^(r) = EV(r). 
2 

The boundary condition for \P(r) is set as 

(0) = 1 and *(oo) = 



. The ground 
(1) 

(2) 



2 



Now a trial function $(r), satisfying the same boundary condition as \I/(r), is introduced 
into the following Schroedinger equation: 

h^V 2 + V (r)]$(r) = Eo$(r). (3) 
The trial function should be chosen such that the difference of the two potentials 

U(r) = V(r) - V (r) (4) 
is small. Defining a Green function D satisfying the following equation[2] 

[~V 2 + V (r)-E ]D = 1, (5) 

then the solution \P(r) can be formally expressed as 

= $-£>([/"- A)tf (6) 
with A = E — i? . Expressing \1/ = $/, then / satisfies 

/ = l-D(U-A)f (7) 

where 

D = <$>- l D<$>. (8) 

Following the same procedure as in [2], in the spherical coordinate system with vari- 
ables (r, 9, (f>) the single trajectory is chosen as the radial coordinate r and correspond- 
ingly we define 



(9) 



hi = 


[(Vr) 2 ] 


-' = 1, 


hi = 


[(W) 2 ] 


~ l =r\ 


hi = 


[(V0) 2 ] 




K = 


hehs = 


r 2 sin 9 



with the volume element 

dr = h^dr d9 d<j>. (10) 

Introducing a step function 

( mi >\ - S for r > r' > , . 

(r 6 r ) = < _ , _ , (11) 

v 1 1 ' \ 0, for < r < r'. v 7 



or 



frieirO = I °' f ° r r>/>0 (12) 

lr|U|rJ -i -1, for 0<r<r'. 
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It can be proved [2] that the Green function is 



= -20$- 2 r- 2 e$V. (13) 

According to the boundary condition of the wave functions we have 

/ = 1 at r = 0, and / = finite at r — > oo. (14) 

Using the definitions of the step function ([11]) and QjjP to express the first and the second 
G-function in (|13"D , respectively, / has the following integral form 

r oo 

/(r) = 1-2 / $- 2 (r')r'- 2 c/r' / $ 2 (r")(t/(r") - A)f(r")r" 2 dr" (15) 



and \l/(r) = $(r)/(r) is the solution of eq.([J). Considering the definitions in eqs.(^) and 
(P70|) , the corresponding energy correction A can be expressed as 

oo 

J § 2 (r)U(r)f(r)r 2 dr 

A = ^ , (16) 

J (& 2 (r)f(r)r 2 dr 
o 

which can be formally expressed as 

A = [Uf]/[fl (17) 

with [F] defined as 

oo 

[F] = J $ 2 (r)F(r)r 2 dr (18) 
o 

2. Iterative solution 

Based on eqs.(^) and ( |I7]) , two iterative sequences {/„} and {A n } are introduced as 

fo = 1, fi, J2, f3, ■ ■ • , fn, ••• 

A = 0, A 1; A 2 , A 3 , A n , ••• (19) 



We require 



/„ = l-D(U-A n )f n -x 

K = [Ufn-l]/[fn-l]. (20) 
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In the explicit form, 

r oo 

f n (r) = 1-2/ $- 2 (ry- 2 dr> J $ 2 (r")(t/(r") - A n ) /„-i(r>" 2 dr" (21) 

r' 

with 

00 00 
A n = J &{r)f n - X {r)U{rydr / J S 2 (r)/n-i(r)r 2 dr. (22) 


For a properly chosen trial function $(r) this iterative procedure will approach the solution 

f(r) = lim f n (r) for all r > 0, 
A = lim A n (23) 

n— >oo 

and *(r) = $(r)/(r), £ = £ + A. 

3. Variational-Iterative Solution of Ground State 
for Yukawa Potential 

The veriational solution of Yukawa potential has been widely discussed in details[3]. In 
ref . [3] using a Hulthen solution as the trial function the result based on variational method 
is compared to the results from different methods. In our paper starting from the Hulthen 
solution we first give the traditional variational result. Using it as the trial function of the 
iterative procedure further improvement could be obtained. The Hamiltonian for Yukawa 
potential is 

H = --^V 2 - g 2 —. (24) 
2m r 

To simplify the notation, based on the following rescaling[3]: 

a 



mg 2 



a 



mg 2 r — > r (25) 
E 



-> E 

mg* 



the Hamiltonian can be expressed as 

1, 



H = --V 2 - 



2 r 

HV(r) = EV(r). (26) 
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Introducing the trial function 



$ A (r) = - — — e - r+1 * Xr (27) 



with a variational parameter A, which satisfies the Schroedinger equation with a Hulthen 
potential: 



l_o Ae 



-Ar 



H x = -^V 2 - 



2 1 - e- Ar 
H x $ x (r) = E A $ A (r), (2* 



^ = -i + iA-iA» (29) 



where 



Based on the variational method the approximate solution of eq.(|26|) could be obtained 
by minimizing the expectation value of H in eq. 



H =- 



J$|(r)(-iV 2 -^)r 2 rfr 



J <&\{r)r 2 dr 
o 

f - 0. ,30) 

Now using the obtained wave function $ A (r) as the zero-th order solution, we introduce 
the formal Green function expression of \&(r) = <& x (r)f(r) in the following way: We can 
define Green function D and -D satisfying 

(H x - E X )D = 1 (31) 
D = ^ 1 D$ X . (32) 

As in Section 2, introducing the energy difference A = E — E x and the potential difference 
of Hamiltonian (|26| ) and (p8| ) 

^) = - — . (33) 

1 — e r 

we obtain expressions for / and A similar to eqs.(0) and (|T7D . The Green function is 

TJ= -20$^ 2 r" 2 e$^r 2 . (34) 



Introducing the iterative series {f n } and {A„} as ( |T9| ) and (|20|), the successive terms can 
be expressed as 
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Ac = 0, / = 1 

Ax = M A = 1 - 27(Z7 - AO 

A 2 = 1^1, f 2 = l-D(U-A 2 )f 1 

[hi 

A n = f n =l-D(U-A n )f n _ 1 . (35) 

L/n-lJ 

Here the trial function for the iterative procedure is the variational wave function <&\. In 
the explicit form, for the first order, we have 

oo oo 

Ax = J <$>\{r)U{r)r 2 dr j J $ 2 x (r)r 2 dr, (36) 

o o 
E x = E x + A t 

and 

r oo 

/i(r) = 1 - 2 J <S>f(r'y- 2 dr' J $\(r"){U{r") - A 1 )r" 2 dr" . (37) 

r' 

It is interesting to compare this first order result to the variational energy and the first 
order perturbation result. From eq.(BD|) 



oo 



H = 2- 



f ${(r)(H x + U))r 2 dr 



J §\{r)r 2 dr 
o 

= E x + A lt (38) 

which is the same as the first order result of the iterative method. If applying the per- 
turbation method, the first order perturbation also gives the same energy. However, to 
further improve the above result the perturbation method would meet problem with the 
unknown excited unperturbed states which are necessary for the calculation of the cross- 
matrix elements in the second order calculation. For the iterative method, in principle, 
there is no problem for the higher order calculation. In general, for the n-th order, we 
have 



A n = J ${(r)f n ^(r)U(rydr/ I ${(r)f n ^(r)r 2 dr. (39) 
o 

E n = E\ + A n , 
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and 

r oo 

f n (r) = 1-2/ <$>-\r'y- 2 dr' J ^ 2 (r")(U(r") - A n )f n ^(r")r" 2 dr". (40) 

r' 

In the next section the numerical results of E\ and Ei are given and compared to the 
result from the variational method for different a of the Yukawa potential in eq.(p6|). 



4. Numerical Results and Discussions 

For Yukawa potential with different parameter a, the variational parameter A of a Hulthen 



potential can be obtained by minimizing the expectation value of Hamiltonian eq.(|26|) 
v.s. A. The values of A corresponding to different a are listed in the first two columns 
of Table 1. For Hulthen potential the upper limit of A to have a bound state is A < 
2. Considering this limitation the highest value of a is taken to be a = 1.15 which 
corresponds to A = 1.9473. For different {a, A} pair the obtained zero-th order energy 
E\, the 1st order energy correction Ai and the corresponding binding energy of the ground 
state Ei = Ex + Aj in the first order iteration are listed in columns 3-5 of Table 1. As 
mentioned in Section 3 the values of E\ for different {a, A} are exactly the same as the 
variational result in ref. [3]. In columns 6-7 of Talbe 1 are given the values of the energy 
correction A 2 and the obtained binding energy of the ground state E 2 = E\ + A 2 after 
the second order iteration. For larger a, the correction in the second order iteration is 
visible. Comparing to the numerical exact solution in the column 8 taken from ref. [4] the 
values after the second order iteration are obviously improved. In principle, the iteration 
could go on to get more accurate result if the integrations in the iteration procedure could 
be performed properly. 

It is possible to solve some excited states based on this method if the proper variational 
solution could be found. For example, the variational solution of the excited 5-states for 
Hulthen potential is known [3], the corresponding solution for Yukawa potential could 
be obtained by the iteration procedure. This method could also be extended to multi- 
dimensional cases as long as one could choose the proper variational solution and find the 
single trajectory to define the Green function and perform the necessary quadratures [1,2]. 
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a 


A 


-E x 


-Ai 


—Ei 


-A 2 


—E2 




0.1 


0.2296 


0.391790 


0.015628 


0.407058 


0.015268 


0.407058 


0.4071 


0.2 


0.4358 


0.305840 


0.020968 


0.326808 


0.020968 


0.326809 


0.3268 


0.25 


0.5327 


0.269121 


0.021796 


0.290918 


0.021798 


0.290920 


0.2909 


0.3 


0.6263 


0.235881 


0.021753 


0.257634 


0.021757 


0.257639 




0.4 


0.8052 


0.178443 


0.019918 


0.198362 


0.019933 


0.198376 




0.5 


0.9750 


0.131328 


0.016756 


0.148084 


0.016789 


0.148117 


0.1481 


0.6 


1.1377 


0.0929452 


0.0131321 


0.106077 


0.0131915 


0.106137 




0.7 


1.2944 


0.0622339 


0.0095107 


0.0717446 


0.0096002 


0.0718341 




0.8 


1.4461 


0.0383507 


0.0062364 


0.0445871 


0.0063520 


0.0447027 




0.9 


1.5936 


0.0206451 


0.00353503 


0.0241802 


0.0036677 


0.0243128 




1.0 


1.7374 


0.00861985 


0.00153843 


0.0101583 


0.00165138 


0.0102712 


0.01029 


1.05 


1.8081 


0.00460320 


0.00083757 


0.00544078 


0.00093748 


0.00554068 




1.1 


1.8780 


0.00186050 


0.00034306 


0.00220356 


0.00041866 


0.00227916 




1.15 


1.9473 


0.000347161 


0.000065800 


0.000412961 


0.000068335 


0.000415496 
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